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ABSTRACT 

A new sample of stars, representative of the solar neighbourhood luminosity func¬ 
tion, is constructed from the Hipparcos catalogue and the Fifth Catalogue of Nearby 
Stars. We have cross-matched to sources in the 2MASS catalogue so that for all 
stars individually determined Near Infrared photometry (NIR) is available on a ho¬ 
mogeneous system (typically K^). The spatial completeness of the sample has been 
carefully determined by statistical methods, and the NIR luminosity function of 
the stars has been derived by direct star counts. We find a local volume luminos¬ 
ity of 0.I2I ± 0.0041/^Qpc“^, corresponding to a volumetric mass-to-light ratio of 
MjLx = 0.31 ± 0.02 AIq/Lj^q, where giants contribute 80 per cent to the light but 
less than 2 per cent to the stellar mass. We derive the surface brightness of the so¬ 
lar cylinder with the help of a vertical disc model. We hnd a surface brightness of 
99L^Qpc“^ with an uncertainty of approximately 10 per cent. This corresponds to 
a mass-to-light ratio for the solar cylinder of M/L^ = 0.34 The mass-to- 

light ratio for the solar cylinder is only 10 per cent larger than the local value despite 
the fact that the local population has a much larger contribution of young stars. It 
turns out that the effective scale heights of the lower main sequence carrying most 
of the mass is similar to that of the giants, which are dominating the NIR light. The 
corresponding colour for the solar cylinder is V — K = 2.89 mag compared to the local 
value oi V — K = 2.46 mag. An extrapolation of the local surface brightness to the 
whole Milky Way yields a total luminosity of Mk = —24.2 mag. The Milky Way falls 
in the range of K band Tuhy-Fisher (TF) relations from the literature. 

Key words: solar neighbourhood. Galaxy: stellar content 


1 INTRODUCTION 

The stellar luminosity function (hereafter LF), i.e. the in¬ 
ventory of stars as a function of their absolute magnitudes, 
is a fundamental property of a stellar population, and has 
wide implications for understanding star formation, and the 
formation and evolution of galaxies. At present, we can de¬ 
termine a complete LF on a star-by-star basis, down to the 
hydrogen burning limit, for stars in the Milky Way only. 
Such LFs have been determined for stars in the solar neigh¬ 
bourhood, and also for open clusters, globular clusters and 
in the Galactic bulge. Conventionally, the LF refers to lumi¬ 
nosities of the stars in the optical bands, such as the V band. 
One of the most influential determinations of the nearby LF 
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is that of Widen (1974), who based the study on the ’Cat¬ 
alogue of Nearby Stars’ ( Gliese|p~969 ), an extensive compi¬ 
lation of data on all stars within 22 pc of the Sun. More 


recently, Flynn et al. (2006) have re-determined the LF us¬ 


ing modern data, including in particular the Hipparcos data 


(Perryman et al. 1997), but finding only small corrections 


to the classic work. Those authors determined not only lo¬ 
cal star number densities (as in Wielenp974 ), but also de¬ 
rived surface densities of stars (i.e. in a column integrated 
through the Galactic disc) as a function of absolute mag¬ 
nitude. In surface density terms, light from the Milky Way 
disc was found to be dominated by emission from main se¬ 
quence stars with My ~ I mag (spectral type A) and by old 
K- and M-giants with My > I mag ( Houk Fesen|I978 ). 


For the analysis of the dynamics and evolution of galax¬ 
ies, it is necessary to determine the masses of the Galactic 
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components. From population synthesis it is obvious that 
the mass-to-light ratio (hereafter MjL^ in solar units) in 
the optical is very sensitive to the contribution of very young 
stellar populations and thus the recent star formation due 


to the dominating light of O and B stars (e.g. Into & Port! 
nari|2013 for the dependence of M/L on the star formation 


timescale in different bands). Additionally, dust attenuation 
reduces the observed light in the optical bands significantly 
at least in late type galaxies and in the Milky Way. As is 
well known, the analysis of galactic rotation curves suffers 
greatly from the uncertainty in M/Ly of the stellar disc. For 
extragalactic systems this is usually done by adopting or fit¬ 
ting a reasonable M/L value of the components. In a new 
approach Martinsson et al. (2013) used dynamic disc masses 


derived by a combination of integral field spectroscopy and 
a statistical scale height determination to derive M/Lk as 
a function of galactocentric radius for a set of galaxies. 

With the increasing number of spatially resolved ob¬ 
servations in the Near Infrared bands (hereafter NIR), NIR 
luminosities are increasingly used for disc mass determina¬ 
tions, primarily in the K band. The variation of M/Lk for 
different stellar populations is much smaller than in the op¬ 
tical bands and the extinction is smaller by a factor of 10 
compared to the V band. Recently near- and far- infrared ob¬ 
servations are used to determine disc masses of extragalactic 
systems (see Martinsson et al. ]M^ | [McGaugh Schomberg| 


2014 and references therein). But there is still no direct 


method to determine for the same galaxy the surface mass 
density and the surface brightness independently. Another 
particular purpose of using M/Lk is to construct models 
of the Milky Way’s structure as constrained by star counts 
in the NIR in various Galactic fields and NIR luminosity 
functions have been measured in many studies, e, 


wood & Jones (1987); Ruelas-Mayorga (1991); Wainscoat et 

al. 

(1992); Lopez-Gorredoira et al. (2002) 

; Picaud, Gabrera- 


Lavers h, Garzon (2003) to probe the Galaxy’s structure. 


Star count methods are more tractable in the NIR because 
of the greatly reduced extinction compared to optical. A lo¬ 
cal K band LF for main sequence stars can also be obtained 
by converting the optically determined LF to the NIR using 
V — K colours averaged over magnitude bins, as in |Mamon| 
& Soneira (1982), but the total brightness in the Kband is 


dominated by giants. 

It is expected that observations in the NIR are better 
suited than those in the optical bands to track the distribu¬ 
tion of stars, especially to probe directly the mass-carrying 
population of late-type main sequence stars (G, K, and M 
dwarfs). Since the solar neighbourhood is the only place 
where we can determine directly both the luminosity and 
the stellar mass density, it is worthwhile to investigate the 
properties of M/Lk based on the best available data for this 
part of the Galaxy. 

In this paper we present a fresh ‘ab-initio’ determina¬ 
tion of the Milky Way LF in the NIR (we use the filter, 
2.2 /xm, throughout the paper but skip the index ’s’ for clar¬ 
ity). Our study is primarily based on Hipparcos stars and at 
the faint end on NIR data to be included in the up-coming 
Fifth Edition of the Gatalogue of Nearby Stars (GNS5; Just 
et al., in preparation). These data have been obtained by 
identifying the GNS5 stars in the Two Micron All Sky Sur¬ 
vey (2MASS |Skrutskie, Gutri Stiening|2006 ). The data set 
drawn from the GNS5 is augmented by samples of stars se¬ 


lected from the Hipparcos catalogue, for which 2MASS data 
is also available. 

Our paper is organized as follows. In section we de¬ 
scribe the construction of our data set and discuss its statis¬ 
tical completeness. In sectionj^we derive the K band LF and 
investigate its implications. We determine the contributions 
by the stars in the various absolute magnitude bins to the 
luminosity and mass budgets of the Milky Way disc in the 
local volume and in the solar cylinder, the principal result 
being that even in the NIR the light is dominated by early 
type main sequence stars and late giants. The consequences 
of this finding and our conclusions are summarized in the 
last section 


2 DATA SAMPLE 


We have constructed our stellar sample by merging two data 
subsets, in order to probe the bright and faint ends of the LF 
respectively. For the bright end of the LF we have extracted 


samples of stars from the revised Hipparcos catalogue (van 
Leeuwen||2007 ), using criteria of absolute magnitude My 
and the parallax limits summarized in columns 1 and 2 of 
Tabled The first three subsamples fall into the survey part 
of the Hipparcos catalogue and are thus volume complete 
by construction. For the last subgroup {d ^ 25pc), we sam¬ 
pled to apparent magnitude V = 10.3 mag, based on the 


determination by JahreiB & Wielen (1997) that the Hip¬ 


parcos catalogue is complete down to absolute magnitude 
Mv = 8.3 mag. The numbers of stars in each subgroup are 
tabulated in the fourth column of Table [T] Stars with rel¬ 
ative parallax errors larger than 15 per cent were excluded 
(except Antares, see below). The numbers of removed stars 
are also given in Table The vast majority of all the sam¬ 
ple stars were then matched to sources in the 2MASS cata¬ 
logue and their absolute Mk magnitudes and V — K colour 
were derived. Only 10 stars did not appear in 2MASS: for 
these stars, K magnitudes were found in the literature or es¬ 
timated from their known spectral types using the relation 


of Koorneef (1983). Since the nominal errors for the bright¬ 


est stars in 2MASS are of the order of 0.2 mag, we have 
compared the 2MASS magnitudes with literature values. We 
find that the differences between 2MASS and literature val¬ 
ues are so small that they can be ignored for our purposes. 
The resulting sample of Hipparcos stars is illustrated as a 
colour-magnitude diagram (GMD) in Fig. It is complete 
to Mk ^ 4 mag with respect to the V band volumes and will 
be used down to the Mk = 3 ± 0.5 mag bin. Note that the 
discontinuities in the GMD reflect Poisson statistics, because 
the subgroups cover significantly different volumes. For the 
analysis we split the stars brighter than Mk — 2.5 mag into 
giants and dwarfs by the dividing line Mk — 4.75(1/ —A) —7 
(see Fig.[^. The resulting sample sizes are also given in Ta¬ 
ble [T] 

The removal of stars with poor parallaxes (i.e. paral¬ 
lax errors > 15 per cent) means that we lose some stars 
which could significantly contribute to the total luminosity, 
in particular at the sparse, bright end of the LF. The bright¬ 
est among the excluded stars are the Ml.5lab supergiant a 
Sco A ‘Antares’ (parallax 5.89 ± 1.00 mas, Mk = —9.93 ± 
0.37 mag), the M3 III giant Hip 12086 (5.16 ± 0.81 mas, 
Mk = —5.72 ± 0.34mag), and the KOIII giant Hip61418 
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Table 1. Volume complete samples of Hipparcos stars. 


My [mag] 

[pc] 

Ntot 

JV(0.15) 


Nd 

^ 0.8 

200 

4560 

104 

2660 

1796 

]0.8, 2.3] 

100 

2039 

21 

541 

1477 

]2.3, 4.3] 

40 

707 

3 

45 

659 

]4.3, 8.3] 

25 

694 

20 

2 

672 


Note. Column 1 lists the absolute V-magnitude ranges, col. 2 the 
distance limits for completeness, col. 3 gives the total number 
of stars, col. 4 the number of excluded stars due to the parallax 
criterion (i.e. the relative parallax error for the stars is greater 
than 15 per cent, except Antares), and col. 5 and 6 the number 
of selected giant and dwarf stars (see also Fig. Q. 


(7.24zb 2.74mas, Mk — —3.28zb0.82mag). We add Antares 
to our sample, because its inclusion corresponds to a factor 
of two in the star number in the Mk = —10 mag bin and it 
adds 7 per cent to the total local luminosity. The contribu¬ 
tion of all other stars is below 2 per cent and is therefore not 
taken into account. The case of Antares also shows, that the 
sampling of the bright end of the LF has large uncertainties 
despite the large volume with 200 pc radius. 

Since we derive distances and also base our absolute 
magnitudes on Hipparcos parallaxes the associated error 
(which we assume to be Gaussian) can affect our sample 
selection. The skewed error distribution after inverting par¬ 
allaxes to distances results in the so called parallax bias 
(Francis 2013) leading to larger mean distances. Another 
effect of the parallax error arising from the volume limit of 
the sample, but of opposite sign, is the Trumpler-Weaver or 
Lutz-Kelker bias ( Trumpler Weaver||1953 Lutz & Kelker 


1973) which increases star counts as the error volume out¬ 


side our limiting distance is larger than inside and more 
stars should scatter in. In the literature these biases are 
often corrected statistically but it was shown that the Lutz- 
Kelker correction depends strongly on the adopted spatial 
distribution of the sample (usually a homogeneous density 
is assumed) and the true absolute magnitude distribution 
(see e.g. [Smith 2003). The latter may be well defined for 
special stellar types like red clump giants or cepheids. In 
our case for the bright Hipparcos stars there is no simple 
way to calculate the Lutz-Kelker correction due to the wide 
spread in luminosities and stellar types. In order to assess 
the impact of these biases on our star counts we select all 
stars in the Hipparcos catalogue. We sample the parallax of 
each star 100 times randomly as a Gaussian with the mean 
of the original star’s observed parallax and the standard de¬ 
viation of its associated error. Then we perform the same 
cuts as before, divide the counts by 100 and compare them 
to our sample. The 40 pc sample decreases by 1 per cent, 
the 100 pc sample increases by 2 per cent, and the 200 pc 
sample decreases by 2 per cent. Relaxing the parallax error 
cuts yield similar corrections. In effect the corrections are 
small and show no clear trend in star counts as well as in 
the luminosity function so that we can safely neglect them. 

A more severe source of error is dust extinction close to 
the Galactic plane. Due to missing 3D maps of the highly in¬ 
homogeneous dust distribution in the solar neighbourhood 
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Figure 1. Colour-magnitude diagram of the Hipparcos stars used 
for this study. The colour coding is as follows - green: stars with 
My ^ 0.8 mag and within 200 pc of the Sun; dark blue: stars 
with 0.8 < My ^ 2.3 mag and within 100 pc; red: stars with 2.3 

< My ^ 4.3 mag and within 40 pc, and light blue: stars with 4.3 

< My ^ 8.3 mag and within 25 pc, respectively, see also table 
In the construction of our LF we have used Hipparcos stars with 
Mk ^ 3.5 mag (dashed line). The solid line cuts our giant sample 
off our dwarf stars. The black cross in the light blue dotted area 
denotes the Sun. 



we cannot correct for this error and only give rough esti¬ 
mates on its magnitude. For that we use the analytic ex¬ 


tinction model from Rybizki & Just (2015) based on Vergely 


et al. (1998) which gives a good estimate for the mean ex¬ 


tinction. It represents a slab of constant dust density up to 
a distance of 55 pc from the midplane with a hole at the 
solar position representing the local bubble. The resulting 
extinction vanishes for stars closer than 70 pc (inside the 
local bubble) and at Galactic latitudes |6| > 52^. The ex¬ 
tinction reaches its maximum of about Ay = 0.19 mag at 
the limiting distance of 200 pc in the Galactic plane. Ap¬ 
plying the same procedure as before but now also correct¬ 
ing for extinction nothing changes for the 40 pc sample. For 
the 100 pc and 200 pc samples the star counts increase by 
2.4 per cent and 6.6 per cent, respectively. Overall, neglect¬ 
ing extinction leads to a slightly underestimated luminosity 
function for bright stars with My < 2.3 mag. This bias is 
stronger for dwarfs at the upper main sequence, because of 
their stronger concentration to the midplane. On the other 
hand the derived luminosity distribution and total bright¬ 
ness in the K band is dominated by giants, which have a 
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Figure 2. Colour-magnitude diagram representing our CNS5 
stars. Shown here are only stars with reliable distance estimates: 
For 3520 stars (black) trigonometric parallaxes are available, 
while for the remaining 365 stars photometric distances were 
adopted (red symbols). In the actual construction of the LF we 
have used CNS5 stars fainter than Mx = 3.5 mag (dashed line). 


larger scale height and therefore being less obscured by ex¬ 
tinction. 

The faint end of the LF has been determined using an 
updated version of the Fourth Catalogue of Nearby stars 
(hereafter CNS5, in prep.; CNS4, Jahreifi Wielen||1997 ). 
All CNS5 stars within 25 pc were cross matched with the 
2MASS catalogue. From the 4622 CNS5 stars within 25 pc 
135 stars got K-magnitudes from other sources. Only 6 close 
binaries were removed at all as well as 9 brown dwarfs and 
one white dwarf below the magnitude limit of the 2MASS 
survey. For the missing stars, iC-magnitudes were obtained 
from the literature, from spectral types, or applying an Mk- 
{y — K) relation based on CNS5 stars with accurate par¬ 
allaxes and reliable colours. The resulting CNS5 sample is 
illustrated in Fig. as a colour-magnitude diagram Mk vs 
{J-K). 


2.1 Completeness and vertical profiles 

Our Hipparcos samples are volume complete by construc¬ 
tion. But beyond about ±40 pc perpendicular to the Galactic 
plane, stellar densities decline with increasing distance from 
the plane, and lead to significant correction factors convert¬ 
ing the mean luminosity density in the observed volume to 


the local luminosity density at the Galactic midplane. It can 
be easily shown that the impact of a (distance dependent) in¬ 
completeness, if present, on the cumulative number of stars 
A^z as function of 2 : is mainly an apparently reduced local 
density no, but the shape Nz/no is essentially unaffected. 
We thus use Nz to determine the vertical density profiles 
and correct for the local volume density at the midplane 
for the ro =100 pc and the 200 pc samples, separately for 
both giant and dwarf stars. Within a sphere of radius ro, 
a constant density is obtained if the cumulative number of 
stars Nz follows the relation Nz = 27rnoz{rQ — jZ) (grey 
dotted lines in Fig. [^. We tested different vertical density 
profiles (linear, exponential with and without a shallow core 
or Gaussian) and investigated the impact of an offset of the 
solar position Z(^ from the midplane. It turns out that the 
profiles of the dwarfs in both samples can be better fit with 


zq = 15 pc. 

For the 200 pc sample, the cumulative star counts devi¬ 
ate significantly from the constant density fit of the inner 80 
and 50 pc for giants and dwarfs, respectively, which yields a 
local number density of 8.8 x 10“^pc“^ for the giants and 
1 per cent less for the dwarfs. Instead the giants are best 
fit by an exponential profile with a local number density 
no = 9.35 X 10“^pc“^ and an exponential scale height of 
456 pc. The dwarfs are better fit by a cored exponential pro¬ 
file with flat density at the midplane and local density of 
no = 9.90 X 10“^pc“^ and an exponential scale height of 
Zexv — 44 pc at large jzj ^ The flat profile near the 
midplane and the small scale height as well as the corre¬ 
sponding half-thickness of /i = 88 pc are expected for B and 
early A stars with velocity dispersions a ~ 5kms“^ (see 
Fig.|9|. The result is also consistent with the vertical den¬ 
sity profile of the A star population derived in |Holmberg 
Flynn| ( [ 2000 ). We note that the midplane densities are 5-10 


per cent larger than the value from the linear fit at small l^j. 
The conversion factors /200 = I (^) from the mean number 
densities (n) = N^n/V 200 (Woo is the volume of the 200 pc 
sphere) to the local number densities no are f 2 oo,g = 1.18 
and f 2 oo,d = 1.85 for giants and dwarfs, respectively. 

Investigating the 100 pc sample yields an almost con¬ 
stant density for the giants. The best fit with an exponential 
density profile yields a local density no = 1.38 x 10“^pc“^ 
and an exponential scale height of 500 pc (see lower panel of 
Fig.§. The dwarfs of this sample are best fit by an exponen¬ 
tial profile with a local number density no = 4.34 x 10“^pc“^ 
and an exponential scale height of 173 pc. The correspond¬ 
ing conversion factors are fioo,g = 1.07 and /ioo,d = 1.23 for 
giants and dwarfs, respectively. 

Our faint star counts are based on the GNS5 and suffer 
from incompleteness, which we assume to depend on dis¬ 
tance only. We have assessed the completeness of the GNS5 
by carrying out radial cumulative star counts in each magni¬ 
tude bin (with a width of 1 mag). These are shown in Fig.[^ 
In a spatial homogeneous sample the cumulative number of 
starsn grows with distanced as n oc d^. In a double loga¬ 
rithmic log(n) vs (—log(7r)) representation of Fig.|^- where 
TT is the stellar parallax - a homogeneous sample would ap¬ 
pear as a straight line with slope 3 (green lines). 

The completeness limit of the GNS5 in each magnitude 
bin can be seen by a deviation of the actual star counts from 
the line of slope 3. From FigJ^we read off the completeness 
limits summarized in Table pPThe magnitude bins Mk = 13 
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Figure 3. Top panel: Cumulative distributions Nz in \z — zq \ of 
the 2 660 giants (red) and 1 796 dwarfs (blue) in the 200 pc-sphere. 
Best fits of the vertical density profiles are shown in green and 
pink lines, respectively. For comparison the grey dotted line shows 
the expectation for a constant density. The differences of model 
and data A^v are also shown in the inset. The slopes at small 2 : 
of giants and dwarfs determining the local number density differ 
by a few per cent only, which is by chance due to the sample 
selection in the 200 pc sphere. The lower panel shows the same 
for the 100 pc sample with 541 giants and 1477 dwarfs. 



log(7r["]) 

Figure 4. Radial cumulative star counts of CNS5 stars in mag¬ 
nitude bins (bin width 1 mag) ranging from = 4 mag and 
brighter to Mx = 10mag (grey points). For clarity log(n) has 
been shifted arbitrarily along the vertical axis for all but the 
Mx = 7 mag sample which is plotted in black. The green lines 
indicate spatially homogeneous samples with the adopted local 
density. 


give within statistical errors consistent results up to an ab¬ 
solute magnitude of Mx — 5 mag. Beyond that magnitude 
the outer shells yield reduced star numbers, because they 
become increasingly incomplete. 


and 14 mag are dominated by white dwarfs. For magnitudes 
fainter than Mx — 14 mag the completeness limit cannot 
be reliably determined. In order to derive a lower limit for 
the star number densities we have considered a volume with 
a radius of 10 pc. White dwarfs in the solar neighbourhood 
are completely sampled out to a distance of d ~ 13 pc from 
the Sun ( Holberg et al.||2008 ). The stellar number densities 
given in the next section have been determined within each 
completeness limit and were then converted to a standard 
(spherical) volume with a radius of 20 pc. 

The luminosity function of the CNS5 stars alone is il¬ 
lustrated in Fig.[^ In order to demonstrate how the volume 
completeness of the CNS5 influences the predictions of the 
LF, we have split up the original volume of the CNS5 (with 
a radius of 25 pc) into spherical shells of 5 pc width. The lu¬ 
minosity functions derived from each shell are over-plotted 
onto each other in Fig. As can be seen here all shells 


3 LUMINOSITY FUNCTION AND ITS 
IMPLICATIONS 

We constructed the local NIR luminosity function 4 >(Mk) in 
terms of star numbers in the 20 pc sphere V 2 Q — 33 510 pc^ as 
described in the previous section by combining the samples 
of stars brighter than Mx — 3.5 mag drawn from the Hip- 
parcos catalogue and fainter stars from the CNS5. We will 
now discuss ^{Mx), the luminosity budget, and the mass 
functions that result. 


3.1 Local luminosity function 

We show the results of our determination of the local lu¬ 
minosity function in Table Given errors indicate the 
(usually dominating) Poisson errors. We first note how 
smoothly the Hipparcos and CNS5 samples join together in 
the Mx — 2 and 3 magnitude bins, see columns 3 and 5. The 
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Figure 5. Illustration of the distance effect on the NIR LF. The 
volume of the CNS5 has been split into spherical shells of 5 pc 
width. The LFs derived from each shell are over-plotted in colour 
coding onto each other. 


Hipparcos sample is volume complete down to Mr = 4mag. 
However, binary stars have been treated more carefully 
in the CNS5 than in the Hipparcos sample as evinced by 
slightly larger star numbers in these magnitude bins (cf. Ta¬ 
ble]^. In column 9 of Table the luminosity function of the 
dwarfs alone (main sequence and turnoff stars) is shown. 
This was derived by excluding giants according to the di¬ 
viding line in the CMD as discussed earlier (cf. Fig. and 
similarly white dwarfs in the bins Mr =11-14 mag were re¬ 
moved. 

The total LF and the contributions of dwarfs and giants 
are shown in the top panel of Fig. Since the number of 
giants is very small, the giant LF is multiplied by a factor of 
100 to make it visible. The dwarf LF shows a clear maximum 
at Mr = 7-8 mag beyond which it drops by nearly an order 
of magnitude at Mr = 13 mag. The faint end of the LF in the 
L, T dwarf regime is unreliable due to the large error bars 
and incompleteness. The shape of log{^{MR)) is roughly 
consistent with the LF observed in the optical bands as is 
illustrated in Fig.[^ bottom panel. Both the NIR and V band 
main-sequence LFs are shown together with the analytical 
ht by Mamon & Soneira (1982) to the V band LF of Widen 


(1974), and its theoretical transformation into the 2.2/xm 


hlter band. As can be seen from Fig. the analytic models 
ht very well in the F - K dwarf regime. At the very bright end 
{Mr ^ 1 mag) the observed LF falls below the theoretical 
prediction. 


3.2 Local luminosity distribution 


In the previous section we have given the luminosity func¬ 
tion 4 >(Mk), i.e. the number of stars in absolute magnitude 
bins. In order to calculate the contribution of each bin to the 
total local luminosity in V 20 , we now multiply ^{Mr) with 
the mean luminosity Lr — of the bin. 


where Mrq denotes the absolute iF-magnitude of the Sun 
Mrq = 3.27 mag ( [Casagrande et a'r]| 2012 ). The result is il¬ 
lustrated in Fig.j^and again the contributions of dwarfs and 
giants are shown. The dwarf luminosity distribution peaks at 
Mr = 2 mag corresponding to F type stars and is analogous 
to the peak at My — 1 mag in the optical luminosity distri- 
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Figure 6. Top panel: Histogram of ^{Mr) for all stars (grey his¬ 
togram), dwarf stars (blue) and giants (red). For clarity ^{Mr) 
of the giants has been enhanced by a factor of 100 in the red 
histogram. The contribution of white dwarfs at Mr 13 mag is 
clearly visible. Bottom panel: Comparison of the main sequence 
LFs in the V band (triangles) and in the NIR (pink circles). The 
blue line is an analytical fit to the optical data and the pink one 
a theoretical transformation of it to the K band, ^ms is given as 
number of stars in the 20pc sphere V20- 


bution ( Flynn et al.|2Q06 ). The contribution of the faint end 
with Mr > 8 mag is completely negligible. The luminosity 
distribution of giants does not show a clear maximum. In 
contrast, the contribution at the very bright end of super¬ 
giants beyond Mr = — 5 mag is dominating in the NIR. Due 
to the low number statistics at Mr < —8 mag the contri¬ 
bution of the brightest end is very uncertain. Red clump 
giants are responsible for the large value at Mr = —2 mag. 
In the V band the luminosity distribution of giants strongly 
declines for My < 1 mag ( Flynn et al.||2Q06 ). 

We hnd for the total K band luminosity pR — 0.1205 zb 
0.0036 LiC 0 pc“^, now converted to the local luminosity den¬ 
sity. The uncertainty is dominated by the bright end of the 
giants. The giants dominate with 0.0971 zb 0.0036 LiC 0 pc“^ 
(80 per cent) and the dwarfs contribute 20 per cent with 
0.0234zb0.0005 Lrqpc~^ . A re-calculation of the V band lu¬ 
minosity of the same sample yields 0.053 Ly(Z)pc~^ (slightly 


smaller than the old value py = 0.056 Lv^ 0 pc of Flynn et 


|aL]|2Q0^ . Combining the local K- and V band luminosity 
results in {V — K) — 2.46 mag in the solar neighbourhood 
(with {V — K)q = 1.56mag). This colour is slightly bluer 
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Table 2. Stellar luminosity function 


Mk 

[mag] 

(1) 

Sp 

[1 

(2) 

^Hip 

n 

(3) 

limits 

[pc] 

(4) 

^CNS 

(*) 

(5) 

err 

(*) 

(6) 

^comb 

n 

(7) 

err 

(*) 

(8) 

'^’ms 

(*) 

(9) 

Mass 

i-^©] 

(10) 

-10 

M3,4 I-II 

0.0024 




0.0024 

0.0017 



-9 

K-M2 I-II 

0.0012 




0.0012 

0.0012 



-8 

M8 III 

0.0145 




0.0145 

0.0090 



-7 

A,G I-II, M6,7 III 

0.0333 




0.0333 

0.0101 



-6 

M4,5 III 

0.0767 




0.0767 

0.0095 



-5 

Ml-3 III 

0.186 


0.5 

0.9 

0.186 

0.015 



-4 

MO III 

0.321 


0. 


0.321 

0.020 

0.0093 

19.6 

-3 

B0,1, K4,5 III 

0.699 


0.5 

0.9 

0.699 

0.030 

0.0389 

15.0 

-2 

K2,3 III 

2.611 


4.1 

1.5 

2.611 

0.098 

0.0981 

10.8 

-1 

B2,3,G8-K0 III 

3.492 


1.0 

0.8 

3.492 

0.157 

0.531 

6.1 

0 

B5-A0, F8-G2 III 

5.512 


7.7 

2.0 

5.512 

0.441 

3.135 

3.46 

1 

A2-A5 

18.49 


16.9 

3.0 

18.49 

1.01 

15.73 

1.98 

2 

F5 

53.95 


48.1 

5.0 

53.95 

2.57 

51.94 

1.47 

3 

F8-G2 V 

112.2 


119.8 

7.9 

112.2 

6.8 

111.7 

1.14 

4 

G5-K4 V 

184.3 

25 

220.2 

10.7 

220.2 

10.7 

220.2 

0.80 

5 

K5-M1 V 

76.80 

25 

261.1 

11.6 

261.1 

11.6 

261.1 

0.63 

6 

M2-3 V 


20 

397.0 

20.0 

397.0 

20.0 

397.0 

0.48 

7 

M4-5 V 


20 

512.0 

22.7 

512.0 

22.7 

512.0 

0.27 

8 



10 

496.0 

63.0 

496.0 

63.0 

496.0 

0.16 

9 



10 

384.0 

55.5 

384.0 

55.5 

384.0 

0.11 

10 



10 

200.0 

40.0 

200.0 

40.0 

200.0 

0.07 

11 



10 

120.0 

31.0 

120.0 

31.0 

112.0 


12 



10 

96.0 

27.8 

96.0 

27.8 

80.0 


13 



10 

144.0 

34.0 

144.0 

34.0 

48.0 


14 



10 

48.0 

19.6 

48.0 

19.6 

16.0 


15 



10 

> 16.0 

11.4 

> 16.0 

11.4 

> 16.0 


16 



10 

> 88.0 

26.6 

> 88.0 

26.6 

> 88.0 


17 



10 

> 24.0 

13.9 

> 24.0 

13.9 

> 24.0 


18 



10 

> 32.0 

16.0 

> 32.0 

16.0 

> 32.0 



Notes: Col 1: absolute magnitude; col 2: spectral type (after [Garwood fc Jones| ( |1987| >); col 3: luminosity function determined from 
the Hipparcos sample; col 4: completeness limits of the CNS5 in abs. magnitude bins 4 to 10 mag, radius for upper limits at fainter 
magnitudes; col 5: luminosity function determined with the CNS5 stars; col 6: errors to col. 5; col 7: luminosity function combined from 
Hipparcos and CNS5; col 8: errors to col. 7; col 9: resulting luminosity function for main sequence stars in the Xband; col 10: main 
sequence star masses. 

(*): The luminosity functions in cols. 3, 5, 7, and 9 as well as corresponding errors in cols. 6 and 8 are given in terms of stars per 20pc 
radius sphere V 20 = 33 510pc^ and per magnitude interval. 


than the value of {V — K) — 2.55 mag determined in model 
A of Just Jahreiss| ( |2Q10| based on a Scalo IMF. 


3.3 Local mass distribution 


We next examine how the NIR luminosity function relates to 
the mass distribution of the stars. For this purpose we have 
multiplied the numbers of stars with mean masses of the 
stars in each Mk magnitude bin. For dwarfs, the adopted 
masses are reproduced in the last column of Table For 
absolute magnitudes 3.1 ^ Mk ^ 9.8 mag we have used di¬ 
rectly the absolute magnitude M^-mass relation by | Henry] 
& McCarthy (1993). For main sequence stars brighter than 


Mk = 3 mag we have used the M^^-mass relation compiled 


by Andersen (1991) which we have transformed to Mk by 


applying V — K main sequence colours. For red giants we 
assume a mass of 1.4 zb 0.1 AJ© as derived by |Stello et al.| 
(2008) from asteroseismic observations. Supergiants may be 


significantly more massive than that ( Schmidt-Kaler||1982 ), 
but are so few that their mass densities are negligible. The 
resulting mass distributions for giants and main sequence 
stars, which we have been able to estimate without recourse 
to population synthesis modeling ( Bell Sz de Jong||2001 Zi- 


betti, Chariot & Rix 2009), are illustrated in Fig. as a 


function of absolute Mk magnitude. The small contribution 
of giants is scaled up by a factor of 10 for visibility. In con¬ 
trast to the luminosity distribution (see Fig. the local 
mass density is dominated by the lower main sequence with 
Mk > 2 mag. 


The local mass density of main sequence and giant stars 
is 0.0315 and 0.00058^4© pc“^, respectively. For the to¬ 
tal stellar mass density we need to add the contribution of 
brown dwarfs and white dwarfs. For brown dwarfs we add 
0.002 Mq pc“^ assuming a 50 per cent incompleteness in 
the observational data of late M, T and L dwarfs. For the lo¬ 
cal mass density of white dwarfs we use 0.0032 zb 0.0003 Ad© 
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Figure 7. Local luminosity distribution in the NIR (grey his¬ 
togram). The contribution by giants is shown with the red line, 
the blue line indicates the contribution by main sequence stars. 



Figure 8. Distribution of the mass density p for the stellar part 
of the Milky Way disc (grey histogram). Stellar mass densities are 
given as a function of absolute Mk magnitudes (blue lines: main 
sequence stars, red lines: giants). For clarity the mass densities 
of the giants have been enhanced by a factor of 10 in the red 
histogram. 


pc“^ ( Holberg et al.|2008 ), which is similar to our finding of 
0.0030 zb 0.0009 Mq pc“^. This way, we find a total local 
mass density oi p — 0.0373^4© pc“^ for the stellar compo¬ 
nent. The total stellar density is slightly smaller tha n in ear- 
her determinations (0.039, 0.044, 0.0415^4© pc“^; JahreiB 


Wielen|1997| [Holmberg fc Flynn|2000[ |Flynn et al.|200^ 

respectively). In all three publications a larger contribution 
of white dwarfs was adopted. Additionally, in [Holmberg 
jFlynnj ( f2000| ) and [Flynn et al. (2006) the mass of the upper 


MS {Mv < 2.5 mag) was overestimated, and in Holmberg & 
Flynn[ ( [2000 ) the contribution of BDs was too high. 


The local mass density implies a mass-to-light ratio 
of the stellar component in the local volume of M/Lk = 
0.31zb0.02 A^©/Lk©. This can be compared with the optical 
mass-to-light ratio of 0.707 M.q/Lvq, which is smaller than 
the value of 0.7b A4 q/Lvq derived in Flynn et al. (2006) 
mainly due to the smaller local mass density. 


Dettbarn, J. Ryhizki 


3.4 NIR Surface brightness of the local disc 

The results presented so far rehect the NIR LF as well as the 
luminosity and mass distributions as functions of Mk for the 
Milky Way in the solar neighbourhood. More representative 
for the entire Milky Way disc is the local surface brightness, 
i.e. the local luminosity distribution multiplied by the ver¬ 
tical scale heights of the stars. The calculations of surface 
brightness and surface density depend on the adopted disc 
model, because for most sub-populations there are no direct 
observations of the vertical density prohle available. It was 


demonstrated in Flynn et al. (2006) that the vertical scale 
heights h of the various stellar populations can be roughly 
estimated by their vertical velocity dispersions cr, because in 
a given gravitational potential the vertical scale heights are 
in lowest order directly proportional to the latter. For an im¬ 
proved estimation of the surface brightness we use a higher 
order ht of the (half-)thickness h{a) dehned by E = 2/ino 
connecting the local volume density no and the surface den¬ 
sity E of the tracer. For a non-exponential vertical prohle of 
the tracer the thickness h differs from the exponential scale 
height. Each stellar subpopulation is approximated by an 
isothermal component in the total potential characterized 
by the local density po = 0.102A4© pc“^ and effective scale 
height zo = 280 pc leading to 

7 / X / N ziaY 

hi^) = V 9 • 

y O ZQ 


with z(a) = 




( 1 ) 


The dotted line in the lower panel of Fig. shows the ap¬ 
proximation for the thin disc. The parameters po and zq 
are chosen to reproduce the values h{a) of the detailed lo- 

(2010) as well as the mass 


cal disc model of Just & Jahreiss 


model used in Holmberg & Flynn 


(2000); [Flynn et al. (2006) 


(triangles in the lower panel of Fig. [9|. The parameters are 
optimized such that the ht function can also be used for the 
thick disc with velocity dispersion a ~ 40kms“^. 

The upper panel of Fig. shows the measured vertical 
velocity dispersions in each Mk bin and Hipparcos group for 
the giants. At the bright end the velocity dispersions scatter 
around a constant value of 17.3 km s“^, whereas the faint 
giants in the 40 pc group show a large scatter and a larger 
mean value of 22.7km s“^. Since the outliers do not con¬ 
tribute much to the total luminosity, we will use the same 
constant ag = 17.3 km s“^ for ah giants. For the upper main 
sequence dwarfs in the Hipparcos samples we show the mean 
values in each Mk bin (weighted by the number densities in 
the 20 pc volume in order to avoid a bias due to the increas¬ 
ing age with increasing V — K). To ah lower main sequence 
stars we have assigned the mean velocity dispersion of G, K 
dwarfs - falling in the Mk = 5 mag bin - from Table 4 of 
Jahreifi & Wielen (1997) because the CNS5 is kinematically 


biased to high proper motion stars at lower magnitudes and 
we expect - and assume here - the same kinematics for ah 
these stars. Their a values are shown as open symbols in 
Fig-i In this figure the blue line shows the analytic fit, us¬ 
ing a shifted error function, of (t{Mk) for the dwarfs, which 
we use to convert the velocity dispersions to the correspond¬ 
ing thickness h. 

Additionally a correction for the thick disc contribu¬ 
tion is necessary. Adopting a standard isothermal old thick 
disc with a velocity dispersion of 40 km s“^ corresponds to 
a thickness h = 1031 pc. In the solar neighbourhood we as- 
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Figure 9. Conversion of velocity dispersion to thickness in 2 
steps. Top: Observed vertical velocity dispersions a for dwarfs 
and giants in each magnitude bin together with the analytic ap¬ 
proximations. For the meaning of open symbols, see text. Bottom: 
Vertical thickness h of the analytic fit used for the transformation 
(blue dotted line), for the best Just-Jahreiss-model (JJ model A, 
solid orange line), the values for the mass model used in [Holm- 
berg Sz Flynn] ( |2000| ) (HF2000, orange triangles) and in [Flynn 
et al.| ( |2006| > (red dashed line and F2006 for giants, red circle). 
For comparison the linear extrapolation based on the total local 
density is also shown (green dashed line). 


sume a thick disc fraction of 10 per cent for all giants and 
for the lower main sequence with Mk >1.5 mag. Due to the 
larger thickness this fraction is enhanced accordingly for the 
surface brightness. The surface brightness of each subpopu¬ 
lation i is given by T^K,i — 2hi Lk.i- 

The resulting surface brightness distribution is shown 
in Fig. |10| where we find similar features as in the lo¬ 
cal luminosity distribution (Fig. 0. The distribution is 
strongly dominated by the bright end of the giants, the red 
clump giants are visible in the Mk — -2 mag bin, and the 
dwarfs peak at the F dwarfs. The integrated surface bright¬ 
ness is Eic = 98.7Licopc“^ = 19.9 mag arcsec“^ in total, 
composed by 16 . 2 Lk 0 Pc“^ (16 per cent) for dwarfs and 
82.5Licopc“^ (84 per cent) for giants. Antares alone has 
added 5.3 demonstrating the uncertainty due to 

Poisson noise for the brightest supergiants. 

A similar determination of the V band surface bright¬ 
ness distribution shows a similar shape, but with a less dom¬ 
inant red giant contribution (cf. discussion above). We find 
29.1 — 22.7 magarcsec”^. This value is 19 per cent 


larger than the value of 24.4 LvqY>^ ^ determined by Flynn 



Figure 10. Surface brightness distribution in the NIR. The blue 
and red lines indicate the contributions by the main sequence 
stars and giants, respectively. 


et al. (2006) arising from inconsistencies in the earlier trans¬ 


formation to surface brightness as can be seen by comparing 
their Figures 2, 5 and 6. In the I band, which [Flynn et al] 
(2006) have used to derive the location of the Milky Way 
with respect to the Tully-Fisher (TF) relation, their Figures 
9 and 10 seem to be correct. Combining the K and V band 
surface brightnesses yields {V — K) = 2.89 mag for the solar 
cylinder. 

The corresponding stellar mass surface density derived 
from the local mass distribution (Fig. [^ including white 
dwarfs and brown dwarfs is 33.3 Al© pc which is slightly 
smaller than the value 35.5^4© pc“^ of jHolmberg Flynn 


(2004); Flynn et al.j ( 2006| ). It implies a Kband mass-to- 
light ratio of M/Lk = 0.34A4©/Lk©. The corresponding 
optical mass-to-light ratio is MjLy — 1.14 A4©/L\/©. 


3.5 Tully-Fisher relation 


The surface brightness of the disc in the solar cylinder can be 
used to estimate the total K band luminosity of the Milky 
Way and compare it with the observed Tully-Fisher (TF) 
relation of extragalactic systems. We proceed similar as in 


Flynn et al. (2006) for the I band. The total disc luminosity 


is approximately independent of the adopted radial scale- 
length Hr in the range of 2.5-5 kpc. For definiteness we 
adopt Hr — 3.5 kpc and find with E^ = 98.7Lic©pc“^ 
at the solar radius of i?o = 8 kpc a total disc luminos¬ 
ity of 8.0 X 1 0^° Lk©. Adding the bulge luminosity of 


1.5x10^ Lkq (Drimmel & Spergel 


2001 


Portail et al. 


2015) 


yields a total luminosity of Mk — —24.2 mag. In Fig. |ll| the 
data point adopting a maximum rotation velocity of the 
Milky Way disc of 220 ± 20kms“^ and an estimated un¬ 
certainty in the total luminosity of 0.2 mag is shown. For 
comparison we plotted two determinations of the TF rela¬ 
tion in the K band based on 2MASS data. IKarachentsev eti 


al. (2002) used edge-on galaxies to derive the isophotal TF 
relation in K band using K 20 = 20 magarcsec”^. We have 
adapted their TF relation from Fig. 8 to a Hubble constant 
of Ho = 70kms“^ Mpc“^. For the Milky Way seen edge- 
on K 20 is at ~ 20 kpc including more than 98 per cent 
of the total light. The Milky Way is 0.2 mag brighter than 
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logfViBM [km s"^]) 


Figure 11. The Milky Way (dot with uncertainty bars, MW) 
with respect to the TF relations of [Karachentsev et al.| ( |2002| ); 
[Masters et al.|p014| >, respectively. Dotted lines show the 1-sigma 
scatter in the TFs. 


the TF relation. For an alternative TF relation we show in 


Fig. 11 also the corrected equations of Masters et al. (2008 


2014) applied to an Sbc type galaxy. In this case the Milky 


Way is 0.2 mag fainter than the TF relation. The systematic 
differences of the two TFs arise mainly on the inclination de¬ 


pendence of the K 20 isophotal luminosity (see also Said et al. 


2015) and the uncertainties in the extrapolation to the total 


luminosity. Deeper K band observations may solve this issue 
in the future. For the Milky Way seen face-on the isophotal 
luminosity would be 0.4 mag fainter than seen edge-on. 


4 CONCLUSIONS 

We have constructed a new sample of stars representative 
for the solar neighbourhood. The data set is based on two 
subsets separated by the intrinsic brightness of the stars. 
The bright part comprises stars drawn from the survey part 
of the Hipparcos catalogue. The faint part Mk > 3.5 mag is 
based on the CNS5, an updated version of the Fourth Cat¬ 
alogue of Nearby Stars. Each star has been identified in the 
2 MASS catalogue so that we have individually determined 
NIR photometry in a homogeneous system for all stars (typ¬ 
ically Ks) available. 

For the Hipparcos stars we have selected volume com¬ 
plete samples in My magnitude bins and corrected for the 
vertical density profile to determine the midplane density. 
These subsamples were then resampled in Mk magnitude 
bins for Mk ^ 3.5 mag. The spatial completeness of the 
CNS5 sample has been carefully examined by cumulative 
radial star counts for each Mk bin resulting in a reason¬ 
able estimate of the local number density as faint as the 
Mk = 14 mag bin, which is well in the brown dwarf regime. 
For fainter stars we derived lower limits for the correspond¬ 
ing star numbers per magnitude bin. All star numbers have 
been converted to a fiducial spherical volume with a radius 
of 20 pc (centered on the Sun) and constant density as a 
measure of the local volume density. For a detailed analysis 
we have separated the giants and white dwarfs from main se¬ 
quence stars including turnoff stars and brown dwarfs. We 
have then determined the NIR luminosity function of the 


stars in the Milky Way disc by direct star counts in our sam¬ 
ple. The K band luminosity function shows a strong maxi¬ 
mum at Mk — 7 — 8 mag. At Mk = 13 — 14 mag the white 
dwarfs dominate the star counts. 

The luminosity function has then been converted to 
the luminosity distribution of the stars by multiplying the 
star numbers with the typical luminosities of the stars in 
each absolute magnitude bin. The resulting (midplane) lu¬ 
minosity distribution is strongly dominated by the very 
bright end of giants and super giants. A secondary peak 
at Mk ~ 2 mag is due to A - F type main sequence stars 
while at Mk = —2 mag the red clump giants stick out. 
At the bright end there is no decline measurable, which 
means that the low number statistics of the brightest su¬ 
pergiants dominate the uncertainty of the total local lumi¬ 
nosity in the K band. We find a total luminosity density of 
Pk — 0.121 zb 0.004 Lk0Pc“^, where the giants dominate 
with a contribution of 80 per cent. Combined with the V 
band luminosity density of pv = 0.053 Lv^ 0 pc“^ we find a 
value oiV — K — 2.46 mag for the colour in the solar neigh¬ 
bourhood. 

We have determined the mass distribution of the stars 
as probed by the NIR luminosity function in the same way. 
Quite contrary to the NIR light, the mass of the Milky Way 
is dominated by K and M main sequence stars. We conclude 
from this discussion that the mass-carrying population of 
stars in galactic discs cannot be observed directly in the 
NIR on a star-by-star basis. The total mass density of the 
stellar component is p = 0.0373 M 0 pc“^, which is about 10 
percent smaller than earlier determinations due to reduced 
contributions by white dwarfs and brown dwarfs. The local 
mass-to-light ratio is then M/L^ = Mq/Lkq^ The 

corresponding corrected value in the optical is MjLy — 
O.71A40/L^0. 


For a comparison to extragalactic systems it is impor¬ 
tant to determine the surface brightness of the disc. We 
have used a detailed vertical disc model to derive the ef¬ 
fective thickness of the stellar populations in the magni¬ 
tude bins and took into account a correction for the thick 
disc contribution with a larger thickness. The resulting sur¬ 
face brightness function shows similar features as the lo¬ 
cal luminosity distribution. The total surface brightness is 
99 LkqPC~^ — 5.4 X 10“® W m“^ with 84 per cent result¬ 
ing from giants. This value can be compared to the K band 
surface brightness of the disc determined from DIRBE data 
after removing all point sources yielding 68 LkqPC~^ (Mel¬ 


chior, Combes & Gould 2007). The difference corresponds 


to the contribution of all supergiants with Mk < —8 mag, 
which seems reasonable. The mass model yields a surface 
density of 33.3 M 0 pc“^ and a mass-to-light ratio for the so¬ 
lar cylinder of M/Lk = PMMq/Lkq- The correspond¬ 
ing optical surface brightness and mass-to-light ratio is 
29.1Lv"0pc“^ and MjLy — IAAMq/L y^^ respectively. 
With the redetermination of the surface brightness we have 
corrected a bug in the earlier determination by [Flynn et al.[ 
(2006), which happened particularly in the V band. An ex¬ 
trapolation of the local surface brightness to the whole Milky 
Way yields a total K band luminosity of Mk — —24.2 mag. 
With standard values of the disc properties the Milky Way 
falls between the Karachentsev et al. (2002) and Masters et 


al. (2014) K band TF relations. 


The stellar population in the solar neighbourhood is 
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strongly dominated by young stars compared to the popu¬ 
lation in the solar cylinder due to the much smaller scale 
heights of the young populations. Nevertheless, the mass- 
to-light ratio in the K band is only 10 per cent larger in the 
solar cylinder. The reason for the luminosity of the present 
day giants - dominating the light in the K band - being a 
rough measure of stellar mass, carried by F, G, and K stars 
of the lower main sequence, is the similarity of their age dis¬ 
tributions. The birth time distribution for the precursors of 
the giants - mainly F and G dwarfs - is spread over the age 
of the disc and similar to that of the F, G, and K dwarfs 
still on the main sequence. As a consequence, the dynamical 
evolution is similar and produces comparable scale heights. 
Thus we conclude that the mass-to-light ratio does not vary 
strongly in disc populations with a long star formation his¬ 
tory and a calibration of the absolute value is provided by 
the solar neighbourhood properties. 

The colours and mass-to-light ratios are consistent with 
the stellar populations derived in the local disc model of 


Just & Jahreiss (2010). Into & Portinari (2013) determined 


mass-to-light ratios and colours for disc populations with 
exponentially declining star formation histories. Our val¬ 
ues are roughly consistent with these models for relatively 
flat star formation histories. An extrapolation from the so¬ 
lar radius to the whole disc would shift the colour and 
mass-to-light ratios slightly dependent on the disc growth 
model. Our K band mass-to-light ratio of the solar cylin¬ 
der of 0.34A4 o/I/kq is very close to the mean value of 


0.31 AIo/Lk© of 30 disc galaxies derived by Martinsson et 

IZIpoBl ). 
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